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^ • Abstract 



We have investigated one and two dimensional (ID and 2D) arrays of tunnel 
junctions in partial Coulomb blockade regime. The absolute accuracy of the 



'^ , Coulomb blockade thermometer is influenced by the external impedance of 

^ , the array, which is not the same in the different topologies of ID and 2D 



arrays. We demonstrate, both by experiment and by theoretical calculations 
in simple geometries, that the ID structures are better in this respect. Yet 
in both ID and 2D, the influence of the environment can be made arbitrarily 
small by making the array sufficiently large. 



Typeset using REVT^ 



Coulomb blockade thermometry (CBT) was invented five years ago [|T| and has since been 
established as accurate and practical means to determine absolute temperature 0. Until 
recently only one dimensional (ID) arrays were discussed. An interesting suggestion to use 
two dimensional (2D) arrays was put forward by Bergsten et al. in [Q, where limitations 
in measurement rate and tolerance to fabrication failures were discussed and shown to be 
even less restrictive than in ID arrays. A large 2D array consisting of 256 x 256 tunnel 
junctions yielded absolute accuracy of better than 0.3 % at temperatures from two to four 
kelvin. This is somewhat better, although statistical variations from sample to sample in 
2D arrays have not been reported, than what has been achieved with ID arrays consisting 
of only 20 junctions in series (0.5 % standard deviation in absolute accuracy). Besides 
considering practical thermometry with sensors having not excessively many junctions, it is 
very interesting to understand the influence of the electromagnetic environment on tunnelling 
in arrays. In the case of a small tunnel junction and weak tunnelling (i?T ^ -Rq, with i?T the 
junction resistance and i?Q = h/Ae"^ ~ 6.5 kQ) this has been treated by the phase correlation 
theory with the harmonic oscillator bath as the environment [^,0, and it has been recently 
extended to describe in detail double junction structures [^0 and long ID arrays [0. A 
general approach which includes strong tunnelling in ID arrays is due to Golubev and Zaikin 
1^ with emphasis on the low temperature limit. Our purpose here is to show experimental 
data on a few topologically different sets of arrays, and to discuss the theory in some of the 
less complicated cases. 

CBT is based on partial blockade of tunnelling in the regime where the charging energy 
of single electrons, Eq = 2^^^^^ (A^ is the number of junctions and C = Ci is capacitance 
of the ith junction in a homogeneous array), the thermal energy, k-oT, and the electrostatic 
energy difference across the array, eV , where V is the bias voltage, compete. The significant 
property is that the full conductance curve, G/Gt^, against V, or in more general terms 
against v = eV/Nk-BT, can be calculated with a universal result for not too large values 
of the ratio u = Ec/k-eT, which is the expansion parameter. Here, Gt = R^^ is the 
conductance of the array at large transport voltages. The basic result is the linear one. 



G/G^ = l-ug{v) (1) 

with known universal corrections for not small values of u. Here g{x) = [a;sinh(a;) — 
4sinh (x/2)]/8sinh (x/2) is the function introduced in Ref. |l[. The main result is that 
the full width at half minimum in this linear regime, V^i/2,0) has the value 

^1/2,0 = 5A39NkBT/e (2) 

with again a known correction AV1/2 = V^i/2 — ^1/2,0 ; which is proportional to the normalised 
depth of the conductance dip, AG/Gt'- 

Ari/2/Vi/2,o = 0.3921 1AG/Gt. (3) 

An important feature of CBT is that it is not significantly influenced by random background 
charges when u < 1. The half width V^i/2, which is proportional to T at small values of u, 
serves as an absolute measure of temperature. AG/G^, in turn, is a secondary thermometric 
parameter, inversely proportional to T (again in the linear regime): 

AG/Gt = Ec/QkBT. (4) 

The formulae (|I| - §) are obeyed quantitatively with high precision in ID arrays consisting 
of a large number of junctions. Yet in simpler structures, with, e.g., A^ < 10 junctions 
in a ID array, the influence of the external impedance on tunnelling is signiflcant 0,^,0, 
and the conductance dip broadens. The origin of this is that Eqs. (|l] - ^ are based on the 
assumption that the arrays are perfectly voltage biased at the ends, which is not the case due 
to the nonzero external on-chip impedance. If we neglect the influence of the environment 
and non-sequential tunnelling, we can, however, analyse our arrays using Eqs. (|l| - §) with 
no extra corrections, if we assume that the structures are fairly uniform. In simple cases, 
like with N junction ID arrays, these results are easy to obtain analytically 0,0]' ^lut with 
2D arrays we have had to resort to Monte Carlo simulations. 

Topologically different ID and 2D arrays which we will discuss here are shown schemat- 
ically in Fig 1. The bias is connected between the left and the right end (busbars). The 



practical A^ junction ID arrays used for Coulomb blockade thermometry consist of M nom- 
inally identical arrays in parallel. The parallel connection does not affect the results of Eqs. 
(|I| - ^), but it lowers the impedance of the sensors by 1/M and this way makes them more 
suitable for practical measurements 0. In the basic "aligned" 2D structure the "equipo- 
tential" islands of the neighbouring chains are connected to each other by a tunnel junction 
whose parameters, i.e., i?T and capacitance C, are equal to those of the rest of the junctions 
in the structure. In the "diagonal" 2D structure we denote the number of connections at 
each busbar by M' . 

As the simplest illustrative example of comparing ID and 2D structures, let us take ID 
and 2D arrays with A^ = 2, M = 2. They are schematically shown in Fig. 1 (d). We denote 
such a ID structure by II and the 2D structure by H. They were fabricated on nitridised 
silicon chips by standard electron beam lithography and two angle shadow evaporation of 
aluminum, with oxidation to create the tunnel barrier in between. The size of the junctions 
was nominally 0.2 x 0.6 yum^, and three different geometries of both types of samples were 
employed to check the influence of the internal structure of the array and of its termination. 
For example, island size was varied from 1 /im up to 10 /im. 

Four II structures and seven H structures were measured at T = 4.25 K. The inset of 
Fig. 2 shows a typical example of the conductance of a II structure and an H structure 
fabricated simultaneously on the same chip. In both cases AG/Gt — 0.01, which implies 
very small corrections from the linear result of Eqs. (|I|) and (^ (see Eq. (^). Although not 
very apparent from the two curves in the inset of Fig. 2, the generic features are: (i) the dip 
of the H structure is wider than that of the II structure, (ii) the shapes of the dips are not 
the same; simple scaling of the two curves does not make them overlap, and (iii) the dips are 
nearly equally high, when the junction size is the same in II and H. As will be discussed in 
what follows, (i) - (iii) are not consistent with the simple sequential tunnelling result with 
zero external impedance. The simple two junction result |]I| with no external impedance 
would predict V^i/2,0 = 4.0 mV at 4.25 K. The main result here, instead, is that V1/2 of the II 
and H structures cluster around different values, both higher than V^i/2,05 i-G-5 at 4.53 ±0.02 

4 



mV (l^i/2/Vi/2,o = 1.13 ± 0.005) for the II and at 4.79 ± 0.05 mV (l^i/2/'i4/2,o = 1-20 ± 0.01) 
for H, respectively. The internal geometry of the array did not influence the results. This 
result already suggests that the corrections to the basic results of Eqs. (0 - ^ are larger for 
2D type structures. 

And, to make a comparison between all different array types of Fig. 1, several samples 
of different types with A^ = 8 were measured. The number of parallel connections was 
M = 9 in ID and in "aligned" 2D; the "diagonal" 2D arrays had M' = 4. Simple theory 
(Eq. (ID) predicts Vi/2,0 = 16.0 mV. Typical conductance curves measured at 4.25 K are 
shown in Fig. 2 for ID and "aligned" 2D arrays. Also for these samples AG/G^ ~ 0.01, 
being well in the linear regime of Eqs. (|l|) and (0). The histogram of V^i/2 of all the 
samples measured is shown in Fig 3 (a). The ID arrays show a +4 % correction to Vi/2,0 
(^1/2 = 16.68 ± 0.04 mV on the average). The "diagonal" 2D arrays have a width of 
Vi/2 = 18.20 ± 0.14 mV, i.e., it is 14 % wider than the theory would predict. The "aligned" 
2D arrays show V1/2 = 19.0 ± 0.27 mV, which exceeds 14/2,0 by 19 %. Thus, at least in 
the case of arrays with a relatively small charging peak {AG/Gt ~ 1 %), the ID structures 
are more suitable for thermometry if judged based on the faster decrease of the size (A^) 
dependent correction observed. Few arrays with smaller junctions were also measured: the 
difference from the simple theory becomes smaller for all arrays, but we have not done this 
systematically enough for quantitative conclusions. Our basic conclusion on A^ dependence 
is further supported by the measurements of A^ = 20, M = 21 and A^ = 30, M = 31 ID and 
"aligned" 2D arrays: the data are collected in Fig. 3. 

To interpret the results we apply the existing phase correlation theory of single electron 
tunnelling ||5|-0. Before doing that we first see the relation of the conductance in structures 
H and II if the external impedance of these structures is assumed to be vanishingly small. 
We can then use the "orthodox" theory of single electron tunnelling |p, and derive the con- 
ductance curve either analytically in the high temperature limit (n < 1) [1|], or numerically 
by a Monte Carlo simulation |TD|,[TT|. The obvious result is, that in both structures the 
conductance curve has the form given by Eq. (||), but the capacitance seen by one of the 



islands to the bias lines is 2C in II, whereas in H it is 8C/3. Thus, Eq of H is smaller by a 
factor 3/4, and so is the depth AG/Gt = u/Q. The width of the dip is not different in the 
two structures H and II. The corrections in V^i/2 observed experimentally are thus due to 
less trivial reasons. 

The most obvious corrections to the basic orthodox result arise either from the higher 
order tunnelling events (co-tunnelling and other non-sequential processes) or from the 
impedance of the environment. The first possibility can be ruled out easily: the higher 
order events are important in the limit of small tunnel resistances (i?T < Rq — 6.5 kQ), 
but our results did not show correlation with the tunnel resistance of the samples, which 
varied from 15 kQ up to 50 kQ. Therefore, we analysed our observations based on the phase 
correlation theory (i.e., perturbatively) in the H and II circuits. 

The tunnelling rates Tj through the junction j in the two directions (±) are given by 

/ + 00 fP 

Here 5Fj is the free energy change in the tunnelling event and Pj{E) = ^^ J^^ dte^'^^^^'^'^'f^^^ 
is the probability density of the tunnelling electron to exchange energy E with its electro- 
magnetic environment, where in turn, Jj{t) is the phase-phase correlation function: 

J,(t) = 2 / L^:^LJl{coth(— -)[cos(^t) - 1] - ^sin(^t)}. (6) 

The essential parameter is the real part of the impedance Ztj seen by junction j. This 
can be obtained as the equivalent impedance of the surrounding circuit in parallel with the 
capacitance of this junction 0. All other junctions are described as pure capacitors, based 
on the assumption of sequential tunnelling. 

To describe the environment of the junction circuit we have used the simple resistive 
impedance. The circuits of II are identical to what has been described in Refs. |^,|^ with 
resistance Re at both ends of the double junction(s), see inset in Fig. 3 (c). Similarly 
we assume that there is impedance Re connected at each termination of the H structure, 
as shown also in the inset. The impedances Re[Ztj(co')] in Eq. (|) can then be easily 



calculated and the result is Re[Ztj(a;)] = -^ i^ij/^ )2 for all the four junctions in II, and 
Re[Ztj{Lj)] = ^ i+(5/4)^(t/L(!)l+('i74)(^M)4 ^°^ ^^^ symmetrically positioned four junctions in 
H, and Re[Ztj(cij)] = x i+(i/4Hcj/<^ y^ ^'^^ ^^^ interconnecting central junction in H. The cut- 
off frequency is defined as Uc = {RcC)~^. We calculated the current voltage characteristics 
of the II and H structures as functions of -Re with two values, C = 6 fF and C = 3 
fF, corresponding to the estimated range of capacitances of the samples measured, and at 
T = 4.25 K. Despite the fact that we used just resistive environment we can draw some, at 
least qualitative conclusions based on the results depicted in Fig. 3 (c), where the half width 
of the corresponding conductance dip from the calculation has been plotted against Re- Since 
the " unintentional" on-chip impedance is typically a fraction of the free space impedance ~ 
377 Q, i.e., on the order of 100 Q, we can assume that we are close to or below the maximum 
of Vi/2 on the Re scale for both II and H structures. Therefore the width of the H structure 
well exceeds that of the II structure. If we assume that we are exactly at the maximum 
in both cases, II and H, even the quantitative agreement between the measured and the 
calculated half widths is good: from the calculation we obtain V1/2 = 4.9 mV and Vl/2 = 4.6 
mV for H and II structures, respectively. The coincidence between experiment and theory 
may be partly accidental, but gives a ground for the observed difference in the measured 
characteristics of H and II. The additional qualitative observation that the influence of 
the environment is most pronounced for very weak Coulomb blockade (AG/Gt ~ 0.01) 
can be understood by noting that the position of the maximum of V1/2 on the Re scale is 
approximately inversely proportional to Eq, but its magnitude is almost unchanged. The 
analysis of the larger 2D arrays in a dissipative environment is markedly more complicated 
and we have not tried that. The larger scatter in the data of 2D arrays in Fig. 3 (b) may 
indicate stronger R^ dependence than in ID, similarly to what is predicted in Fig. 3 (c) for 
N = 2. 

In summary, we have investigated the influence of the cross connections in arrays of small 
tunnel junctions on their current voltage characteristics. We have observed that the multiply 



linked circuit structure in 2D arrays leads to strong influence of the external electromagnetic 
environment on the tunnelling junction imbedded in the array. This affects drastically 
the thermo metric properties of not very extended junction arrays in Coulomb blockade 
thermometry in favour of ID structures. 

We thank Antti Manninen, Andrei Zaikin, Per Delsing, Jari Kinaret, Klavs Hansen and 
Tobias Bergsten for discussions, and the National Graduate School in Materials Physics for 
support. 



REFERENCES 

[1] J.P. Pekola, K.P. Hirvi, J.P. Kauppinen, and M.A. Paalanen, Phys. Rev. Lett. 73, 2903 
(1994). 

[2] K.P. Hirvi, J.P. Kauppinen, A.N. Korotkov, M.A. Paalanen, and J.P. Pekola, Appl. 
Phys. Lett. 67, 2096 (1995); J.P. Kauppinen and J.P. Pekola, Phys. Rev. Lett. 77, 
3889 (1996); Sh. Farhangfar, K.P. Hirvi, J.P. Kauppinen, J.P. Pekola, J.J. Toppari, 
D.V. Averin, and A.N. Korotkov, J. Low Temp. Phys. 108, 191 (1997); J.P. Pekola, 
J.J. Toppari, J.P. Kauppinen, K.M. Kinnunen, A.J. Manninen, and A.G.M. Jansen, J. 
Appl. Phys. 83, 5582 (1998); J.P. Kauppinen, K.T. Loberg, A.J. Manninen, J.P. Pekola, 
and R.V. Voutilainen, Rev. Sci. lustrum. 69, 4166 (1998). 

[3] T. Bergsten, T. Claeson and P. Delsing, J. Appl. Phys. (1999); T. Bergsten, T. Claeson 
and P. Delsing, Proceedings of the 22nd international conference on low temperature 
physics, LT22, to be pubhshed (2000). 

[4] M.H. Devoret, D. Esteve, H. Grabert, G.-L. Ingold, H. Pothier, and C. Urbina, Phys. 
Rev. Lett. 64, 1824 (1990); Yu. V. Nazarov, Sov. Phys. JETP 68, 561 (1989) [Zh. Eksp. 
Teor. Fiz. 95, 975 (1989)]. 

[5] G.-L. Ingold and Yu.V. Nazarov, in Single Charge Tunneling, NATO ASI Ser. B, Vol. 
294, edited by H. Grabert and M.H. Devoret (Plenum, New York, 1992), Ch. 2. 

[6] P. Joyez and D. Esteve, Phys. Rev. B 56, 1848 (1997) [E 58, 15912 (1998)]. 

[7] Sh. Farhangfar, A.J. Manninen, and J.P. Pekola, Europhys. Lett. 49, 237 (2000). 

[8] D.S. Golubev and A.D. Zaikin, Phys. Lett. A 169, 475 (1992). 

[9] D.V. Averin and K. K. Likharev, in Mesoscopic Phenomena in Solids, edited by B. L. 
Altshuler, P. A. Lee, and R. A. Webb (Elsevier, Amsterdam, 1991), Ch. 6. 

[10] N.S. Bakhvalov, G.S. Kazacha, K.K. Likharev, and S.I. Serdyukova, Sov. Phys. JETP 



68, 581 (1989) [Zh. Eksp. Teor. Fiz. 95, 1010 (1989)]. 
[11] K.P. Hirvi, M.A. Paalanen, and J. P. Pekola, J. Appl. Phys. 80, 256 (1996). 



10 



FIGURE CAPTIONS 

Fig. 1. Schematic presentation of the different array structures investigated: (a) ID, 
(b) "ahgned" 2D, (c) "diagonal" 2D arrays, and (d) the double junction II (left) and the 
coupled double junction H (right) structures. 

Fig. 2. Measured normalised conductance G/Gt vs. bias voltage V for the topologies 
(a) and (b) of Fig. 1 with N = 8 and M = 9, measured at 4.25 K. Inset shows the 
corresponding conductance curves of the II and H structures in Fig. 1 (d). 

Fig. 3. Summary of the experimental and theoretical results: (a) Histogram of the 
measured V1/2 of the samples with N = 8, M = 9, and M' = 4 (at 4.25 K). The filled black 
bars are for ID arrays, the grey ones are for "diagonal" 2D arrays, and the open ones for 
the "aligned" 2D arrays, (b) the half width of the conductance curve for ID (filled circles) 
and "aligned" 2D (triangles) arrays against A^. The "crossed" 2D data (A^ = 8, M' = 4) 
are shown by a cross. The lines are just guiding the eye. In (c) we show the circuit models 
used in the calculations, and the theoretically obtained half widths for the two structures 
II and H as functions of the external impedance (resistance) Re- In all cases T = 4.25 K, 
and all the junctions are identical with C = 6 fF (C = 3 fF) for the solid (dashed) lines. 
Changing either G or T simply scales the values on the horizontal axis; the main parameter 
in the calculation is [h/{ReG)]/{kBT). 
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